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The properties of the observed Higgs boson with mass around 125 GeV are constrained by a 
wealth of experimental results targeting different combinations for the production and decay of 
a Higgs boson. In order to assess the compatibility of a non-Standard Model-like Higgs boson 
with all available results, we present Lilith, a new public tool that makes use of signal strength 
measurements performed at the LHC and the Tevatron. 


I. INTRODUCTION 

The discovery of a new particle with properties compatible with those of the Standard Model (SM) Higgs 
boson and mass around 125 GeV at CERN’s Large Hadron Collider (LHC) [TJ |2] was a major breakthrough. 
Detailed information on the properties of this Higgs boson was obtained from the measurements performed 
during Run I of the LHC at 7-8 TeV center-of-mass energy E3I1], making it possible to test a large variety 
of models of new physics. The results of the Higgs searches at the LHC are usually given in terms of signal 
strengths, p, which scale the number of signal events expected for the SM Higgs boson. The signal strength 
measurements can be combined and used to constrain new physics models in which i) the signal is a sum of the 
processes that exist for the SM Higgs boson, and ii) the Higgs couplings are defined through a simple rescaling 
of the corresponding SM ones. In this talk we present a new public tool for performing a global fit to all 
available signal strength measurements, Lilith [5]. It is a library written in Python that can easily be used in 
any Python script as well as in C and C++/R00T codes, and for which we also provide a command-line interface. 
As a primary input, it uses signal strengths results for which the fundamental production and decay modes 
have been unfolded from experimental categories. The experimental results are stored in XML files, making it 
easy to modify and extend. The user input can be given in terms of reduced couplings or signal strengths, and 
is also specified in an XML format. 


II. LIKELIHOODS AND PARAMETRIZATION OF NEW PHYSICS 

Each event category (corresponding to a given set of selection criteria) is sensitive to different combinations 
for production and decay of the SM Higgs boson, with various efficiency factors. Given this information, signal 
strength measurements can be combined and re-expressed after unfolding of the fundamental production and 
decay modes as 


p(X,Y) 


a(X)B(H -> Y) 
a SM (X)B SM (H ->■ Y) ’ 


( 1 ) 


for the different production modes X £ (ggH, VBF, VH = WH + ZH, ttH) and decay modes Y £ ( 77 , ZZ *, 
WW*, bb, tt , ...) of the SM Higgs boson. Results have been systematically given in terms of p(X, Y) for the 
main decay modes probed by the ATLAS and CMS collaborations. In particular, results were often published 
in the plane (^i(ggH + ttH, Y), /r(VBF + VH, Y)), see Fig. [l] 

This information is taken into account in Lilith in the following way: for two (combination of) production 
and decay processes (X,Y) and (X', Y'), in the Gaussian approximation the experimental likelihood reads 

-2 log L(p) = (p-p,) T C~\p-p), (2) 


Wher<i ^ = « XDMX-.Y-)) and C- = ^ J is the inverse of the covariance matrix. Tire list of 

experimental results present in the database of Lilith can be found in [5] (see also [6] for updates). It includes 
the results from direct searches for invisible decays of the Higgs boson through VBF and ZH production. 
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FIG. 1: Two dimensional ATLAS (left) and CMS (right) results in which the fundamental production modes are unfolded 
from experimental categories, in the plane (/r(ggH + ttH, Y), /x(VBF + VH, Y)) for V £ ( 77 , ZZ* , WW*, tt) [ 3 ( H ]. 


Whenever available, full 1- or 2-dimensional likelihood functions are used as they make it possible to go beyond 
the Gaussian approximation. All results are combined into a global likelihood function L defined as the product 
of the individual likelihood functions. 

Results given in terms of signal strengths, as in Eq. ([I]), can be matched to new physics scenarios with the 
introduction of factors Cx and Cy that scale the amplitudes for the production and decay of the SM Higgs 
boson, respectively, as follows: 


/r(X, T ) ^ Y ) ' (3) 

The factors Cx and Cy can be identified to (or derived from) reduced couplings appearing in an effective 
Lagrangian Hi- In particular, the factors CggH and C 77 —scaling gluon fusion and the decay into two photons, 
respectively—can be computed from the reduced couplings at LO or NLO QCD accuracy in Lilith. It is also 
possible to specify decays into invisible and/or undetected particles. Given this information, Lilith computes 
the signal strengths (alternatively they can be given directly as input) and evaluates the global likelihood 
function. Contraints can then be put on new physics scenarios in a frequentist or Bayesian approach. 


III. QUICK INTRODUCTION TO LILITH 

Lilith is a library written in Python for constraining models of new physics against the LHC results. The 
code is distributed under the GNU General Public License v3.0. All necessary information on how to download 
and install it can be found on its official website |Bj. Lilith requires Python 2.6 or more recent, but not the 
3.X series. The standard Python scientific libraries, SciPy and NumPy, should furthermore be installed. 

Lilith can be used in three different ways: 

1. from a a Python code or an interactive session of Python, using the Lilith application programming 
interface (API); 

2. calling the command-line interface (CLI) run_lilith.py from a shell; 

3. or using the Lilith interface to C and C++/R00T. 

All three methods are documented in 0. A minimal example of use of the API is now presented: 

from lilith import * 
leal = LilithO 
leal.readexpinput() 

leal.readuserinputfile(’userinput/example_mu.xml’) 

leal.computelikelihoodO 

print ; -21og(likelihood) =’, lcal.l 
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FIG. 2: Fits to the benchmark scenario ( Cv,Cf ) using data from the ATLAS (left) and CMS (right) collaborations [3j[4]. 
The red and orange filled surfaces correspond to the 68% and 95% CL regions obtained by the ATLAS and CMS 
collaborations while the corresponding dashed lines show the Lilith results. The black star indicates the position of 
the Lilith best-fit point, the blue diamond is the ATLAS or CMS best-fit point and the white circle shows the SM 
prediction. 



After importing the Lilith library, an object leal is instantiated. The list of experimental results as well as the 
user input in XML format are then read with the methods readexpinput and readuserinputf ile, respectively. 
Finally, the likelihood is evaluated with the method computelikelihood and the result is printed on the screen. 

The XML user input hie consists in either reduced couplings or signal strengths (possibly with additional 
invisible or undetected decays) for one or several states contributing to the signal. The root tag should be 
<lilithinput>, and is followed by either <reducedcouplings> or <signalstrengths>. For instance, scaling 
factors associated with the H —»• bb and H —> 77 processes can be specified as follows, 

<C to="bb">l.2</C> 

<C to="gammagamma">0.8</C> 

The complete list of available attributes can be found in the manual mm Alternatively, the signal strengths 
can be provided as input. For instance, taking /r(ggH, WW *) and bb) it would read 

<mu prod="ggH" decay="WW">l.3</mu> 

<mu prod="ttH" decay="bb">l.l</mu> 

Finally, with reduced couplings a <precision> tag can be given, with possible values being "L0" or "BEST-QCD" 
for the computation of loop-induced couplings at LO or NLO QCD accuracy, respectively. 

We also provide a CLI as well as an interface to C and C++/R00T in order to easily use Lilith outside of a 
Python script. All details and documentation are given in the manual [5]. 


IV. VALIDATION AND EXAMPLE 

The Higgs likelihood used in Lilith should be validated against official LHC results as approximations are 
made when reconstructing the full likelihood from 1- or 2-dimensional results. To this aim, we reproduce global 
fits to reduced couplings in various scenarios considered by ATLAS and CMS. An example is given in Fig. [2] 
It shows very good agreement with the results from both experimental collaborations. More examples can be 
found in the manual [5j. 

Several examples are shipped with Lilith (for the complete list, see [5]). In the Python code 
CVCF_ldprofile .py, 1-dimensional constraints on Cy and Cp are derived from a global fit to ATLAS and 
CMS results in the (Cy, Cf) benchmark scenario. Results are shown in Fig. [ 3 ] This example uses the library 
iminuit [2], a Python implementation of the MINUIT nni library, in order to minimize —2 log L and derive the 
profiled likelihood around the minimum. The plotting library matplotlib m is used to produce the figures. 
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FIG. 3: 1-dimensional likelihood profiles of Cv (left) and CV (right) from a global fit of the benchmark scenario (CV, Cf). 


V. CONCLUSIONS 

Using all available information from the ATLAS and CMS collaborations to construct a likelihood is a non¬ 
trivial task. To this aim, we provide a new public tool, Lilith. Lilith is a library written in Python, and 
for which we provide an API as well as a command-line interface and a basic interface to C and C++/R00T. The 
experimental results are read from a database in XML format that is shipped with the code and which is easy 
to modify and extend. Lilith uses as a primary input results in which the fundamental production and decay 
modes are unfolded from experimental categories. 

New physics can be parametrized in terms of reduced couplings, or signal strengths directly, which are given 
as input to Lilith in XML format. If needed, scaling factors for the loop-induced processes and VBF production 
are computed taking into account QCD corrections. The likelihood is evaluated from a set of experimental 
results and given as output. The Higgs likelihood of Lilith obtained from the latest measurements at the LHC 
has been validated against ATLAS and CMS results and can be used to constrain new physics. 
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